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Abstract 

The motion of spectator quarks in decay of a beauty hadron is a nonperturbative effect which 
can usually be neglected. We find that the motion in some decay channels, which contribute total 
decay widths of beauty hadrons, can not be neglected. The contributions from these decay channels 
to decay widths are proportional to certain averages of the squared inverse of the momentum carried 
by a spectator quark. This fact results in that these contributions, suppressed by l/m^ formally, 
are effectively suppressed by l/m&. We find these contributions can be factorized into products of 
perturbative coefficients and nonperturbative parameters. We calculate these coefficients and define 
these nonperturbative parameters in terms of HQET matrix elements. Since these parameters are 
unknown, we are unable to give numerical predictions in detail. But with a simple model it can be 
shown that these contributions can be large. 
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Although the heavy quark effective theory(HQET) is very successful to predict various properties of 
hadrons containing one heavy quarkpQ, but it is still difficult to predict the experimentally observed ratio 
of decay widths of S-meson and 6-navored baryon. The experimental results for lifetime ratios of beauty 
hadrons are|2j: 

T ^ B+ ] = 1.074 ± 0.014, 4^4 = 0.948 ± 0.038, = 0.796 ± 0.052. (1) 

r(B d ) r{B d ) r(B d ) V ; 

In HQET a systematic expansion in l/m& is employed to give theoretical predictions. At the leading 
order, the decay width of a beauty hadron equals the decay width of a free 6-quark, hence the above 
ratio should be one, if higher-order correction is neglected. It is clearly that from Eq.(l) the ratio with 
Aft deviates from one substantially. 

rder result starts at order of rrv b 
hadron, denoted as Hb, takes in general the form 



The correction to the leading-order result starts at order of m b 2 . The decay width of a b- flavored 



r ( ft) = r (6 ) 1 + Y><^!» , (2 ) 

V n=2 b / 

where c n is a perturbative coefficient, O n 's are operators of HQET, whose matrix element represents 
nonperturbative effects in the decay. So far only contributions from local operators are considered. It 
should be noted that the correction at order of m b does not exist. It is expected that spectator effects 
will results in the difference between life times of 6-flavored hadrons. Spectator effects appear at order 
of m^ 3 . These effects have been studied in (HIE] at tree-level of QCD. Next-to-leading order corrections 
have been studied [El El Cj. In these studies the nonperturbative effects at order of m^ -3 are parameterized 
with matrix elements of four quark operators. These matrix elements have been studied with lattice 
QCD 8 or with sum rule techniques(e.g., see [H]). With obtained values of matrix elements and effects 
at next-to- leading order the ratio becomes [7j: 

t(B + ) t(B s ) r(A{,) 

-—i = 1.06 ± 0.02, -4ttt = 1-00 ± 0.01, -4ttt = 0.90 ± 0.05. (3 

T{B d ) T{B d ) T{B d ) 

This prediction is closer to experiment, but there is still a discrepancy at order of 10% for A&. Other 
attempts to explain the ratio in Eq.(l) can be found in ^01 • It is interesting to note that the effects at 
order of m b 5 studied before consist of local operators of four quark fields, in which two fields are fields of 
HQET for 6-quark, while other two are for light quarks. These light quarks can be either as a part of Hf> 
as a bound state, or they generate nonperturbatively soft decay products. Since the operators are local, 
the light quarks represented by the two quark fields carry zero momenta in the decay, i.e., the motion of 
spectator quarks is neglected. 

The formula in Eq.(2) is based on a factorization in HQET, in which nonperturbative effects are 
factorized from perturbative effects. In general, one can expect that similar factorization formula can be 
obtained for inclusive productions of H&, where nonperturbative effects related to Hb are parameterized 
by matrix elements of local operators, in corresponding to those matrix elements of {Hb\O n \Hb) ■ At the 
leading order of m^ 1 the production rate of Hb can be factorized as a product of the production rate of a 
free 6-quark with a matrix element of HQET. The matrix element can be interpreted as the probability 
for the transition of the 6-quark into Hb- Inclusive productions based on such a factorization have been 
studied, predictions at leading order of m b l have been made for production at e + e~ colliders and for 
polarization of heavy vector meson [TTlll2| , In these cases a good agreement with experiment were found. 
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Recently, such a factorization was employed to explain the asymmetry between production rates of 
D + and of D~ in their inclusive productions, and also asymmetries for other heavy flavored hadrons|13| 
1141 1151 ITS] . In these works contributions from quark recombination to production rates are studied, in 
which a heavy flavored hadron like Hb is produced by combining a 6-quark with a light antiquark q. 
Including these contributions the asymmetry can be explained T3~] IT U 115) . The nonperturbative effect 
of the recombination can be represented as matrix elements of four quark fields in the production rate. 
But it is found that the momentum of q can not be taken as zero, because it will generate a type of 
infrared singularities in the production amplitude of bq. This type of singularities can be regularized by 
noting the fact that the light antiquark inside Hb carries a small momentum at order of Aqcd and its 
effect is nonperturbative. One needs new matrix elements beside these matrix elements of local operators, 
corresponding to those in Eq.(2), to incorporate this nonperturbative effect. These new matrix elements 
are found to be nonlocal matrix elements, whose definitions are given in |16j . It turns out that the 
contributions to the production rates due to quark recombination is proportional some averages of the 
inverse of momenta carried by the light antiquark. This results in that the contributions are significant 
and the large asymmetry observed in experiment can be explained in this way. The effect in quark 
recombination in production c-flavored barvon[T7] and in production 6-flavored jet in Z°-decay^S] has 
also been studied. 

In this letter we point out similar contributions also exist as corrections to lifetime ratios. These 
contributions formally are suppressed by m^ 3 , but they are proportional to the square of the inverse of 
the momentum carried by a light antiquark in Hb or by a light antiquark which generates soft decay 
products nonperturbatively. Therefore these corrections are only suppressed by m^ 1 effectively. Similar 
situation also appears in the decay B — * jiv, where the decay amplitude is proportional to a certain 
average of the inverse of the momentum carried by the light quark inside the B meson. This decay mode 
has been studied in detail in |19| I2flj. In this letter we will study this type of contributions to lifetimes 
of beauty hadrons. 

The effective weak Hamiltonian for the decay of Hb at the scale [i = nib is: 

G 

H e ff = ~J=VcbJ2 { c l( m bWugQLYu L C L J^b L + V* q q L YcLCL7^b L } 
* q=d,s 

+c 2 (m 6 ) [Vu q CLl^u L q L ^^b L + V* q c L ^c L q L ^^b L } } , (4) 

where we neglected the suppressed transition b — ► u and qi = 5(1—75)5". c\ and c 2 are Wilson coefficients. 
The above mentioned contributions come from these decays at the tree-level: 

H b (P) - c ( Pl ) + q( P2 ) + G(k)+X, 
H b {P) -» c { Pl )+u{p 2 )+G(k)+X, 

H b (P) - c( Pl ) + c(pa) + G(k) + X. (5) 

where momenta are given in brackets. We will not consider a c-quark as a spectator, hence the decay into 
a lepton pair with a gluon and other unobserved states X will not lead to contributions we consider. In 
general the gluon can have any possible momentum. If the momentum is large, one can use perturbative 
QCD. If the momentum is small, the contributions can have an infrared singularity, reflecting the fact 
that the gluon can not be taken as a perturbative gluon. However, it turns out that the contributions 
we are interesting in are free of infrared singularities. This will be discussed in detail. We will give some 
detail for calculation the contribution from the process Hb — > c + q + G + X . 
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With the effective Hamiltonian the decay amplitude can be written as: 



T = J 72^M?i>*0 / d^e^HXlu^x^mHb), (6) 

where i,j stand for spin- and color indices. is the scattering amplitude for u(qi) + 6(92) — ► c(j>i) + 
q(P2) + G(k) in which the quarks in the initial state are off-shell in general and q\ + q2 = Pi + P2 + k. 
With the translational symmetry the contribution to the decay width can be written: 

■(-1) J d l x 1 d l x 3 d t x l ^' '•-^ I '-""'(H b lh(x i )u i (x 1 )u t (x 3 )b j (0)lH l }, (7) 

with q4 = pi + P2 + k — q%. In the above we exchanged the order of u-quark fields and this gives the — 
sign in the second line. The integration measure for the phase space of three particles cqg is denoted as 
dT cqg . We use nonrelativistic normalization for the state and 6-quark. The above contribution can 
be illustrated by Fig.l., where one of 16 diagrams is shown explicitly. 




Fig.l 

Figure 1: Diagrams for the contributions to the decay width. Other diagrams are obtained by changing 
attachments of the gluon. The thick line is for the 6-quark. The broken line is the cut, the black box 
represents nonperturbative effect in the decay, its expression is given in Eq.(9). 

For 6-quark field b(x) one can use the m^ 1 expansion: 

b( x ) = e - im » v - x (h(x) + •••), b{x) = e +im » v - x (h{x) + •••), (8) 

where v the four velocity of Hi,. The • • • represent higher orders in m^ 1 , which can be neglected in this 
letter. With v any vector can be decomposed as B^ = v ■ Bv^ + B^_ with B± ■ v = 0. Taking the 
leading term the Fourier transformed matrix element becomes 

J d A Xl d A x z d A x± e^- 1 ~^- :E3 -^^(^|6K^4)^(xi)n fc (x3)6 J (0)|if 6 ) 

= J d 4 Xl d 4 x 3 d 4 Xi e^ 1 - Xl - i93 - X3 -^ 4 -( 94 - m ^)(i7 6 |^(x4)u i (xi)u fe (x 3 )^(0)|i7 6 ) + • • • . (9) 
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Since we extract the large momentum m^v by using the expansion in Eq.(8), the space-time dependence of 
the matrix element in the second line in Eq.(lO) is controlled by the soft scale Aqcd- The x-dependence 
of h fields can be safely neglected. If we can neglect the x-dependence of the light quark fields, then the 
contribution will proportional to matrix elements of local four-quark operators, which appear at order 
of m^ 3 in Eq.(2). This implies that the light quark q will carry zero momentum. But the amplitude A 
with the zero momentum is divergent or proportional to the inverse of the light quark mass m u if we 
do not neglect this mass. This divergence represents some new nonperturbative effects which can not 
be represented by local four quark matrix elements. To study the divergence, we introduce a light-cone 
coordinate system in which the two light-cone vectors are I and n respectively and I ■ n = 1. In the 
light-cone coordinate system the emitted gluon has the momentum fc M = \/2kPl^. The expansion in q\ 
reads: 

A ij (q 1 ,m b v) = -Tij H , (10) 

qi • / 

where • • • stand for higher order terms of q\ and we set qi = m^v. Tj reads: 

Ta = i^^{ci[«(pi)7^ 

-c 2 [u( P2 )r(l - 75)] j [u(Pi)T°7, t (l - 75)7 • e*(k)-y • 
= [T'l-% (11) 

where e*(k) is the polarization vector of the gluon and we used k ■ s*(k) = 0. Tj only gets contribution 
from Fig.l., it does not depend on q±. It should be noted that the matrix T can always be written in the 
form as in the last line of Eq.(10) for our processes in Eq.(5). Keep only the leading term for Aij, some 
integrations in Eq.(7) can be performed. We obtain: 

= / ^„ J 'hi'^m^, — pi — P2 — ( 7 «tv) h 

■(-ml) J du^duz e*"""">-™™""{if i |/.,(0)[7"«]i(w()[»7"]*("3()'>j(0)|« t ). (12) 

where we have used qi-l = rjirrib for i = 1, 3 and moved 7" = j-l. into the matrix element. It is interesting 
to note that the Tij is finite when the gluon carries null momentum. This indicates that the decay width 
will be free from infrared singularities when we perform the phase-space integration. If one keeps the 
next-to-leading order in q\ or q%, infrared singularities will appear, but these singularities may be cancelled 
by some virtual corrections partly and absorbed into four quark matrix elements. The contribution at 
this order will be suppressed by (M^ — mi) /mi in comparison with that from the leading order. In our 
approximation the decay width will be proportional to the integral of the Fourier transformed matrix 
element. The Fourier transformed matrix element can be parameterized as: 

1 1 ihn lh, ' : 1 ' U:i ' U ' 3 e^^-^^^i^KoJb-^iC^iO^-lfe^O^WI^) 



m b J Vi V3 J 2tt 2ir 

= \[(Pv) jk (Pv) u 4 U ; 1} - {isPv) jk (Pvish vtf - (tfp v ) jk (p v jt,) u 

- (rtl5Pv) jk (PvWhA^} + \[{P v T*) jk (P V T«) U S^ _ {l5 P vT «) jk (P V ^) U V^ 8) 

~ (l" T PvT a ) jk {PvlT»T a ) u V£ 8) - {l^T a P v ) jk {P vlW T«) u A^f], (13) 
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with 



7t 



7^ — v ■ JV' 



1 + 7 ■ v 



(u,l) 



(14) 



in 



The above matrix elements are defined in the rest frame of H b . The eight coefficients S Hi 

Eq.(13) are dimensionless, their values depend on H b . These coefficients multiplied with the factor V2v° 
are Lorentz covariant because we take nonrelativistic normalization for the state. We work in the light- 
cone gauge I ■ G = 0. In other gauges gauge links should be added into the above matrix element between 
quark fields to make it gauge invariant. If we replace in the integral the factor (ryi^) -1 with a constant, 
the eight coefficients will be at order of K 3 q CD /m 3 . Taking this fact into account, the coefficients are at 
order of m 2 /(MH b — m b ) 2 ■ Aq^/to 3 , ~ A-QCD/mb- Therefore, the contribution to the total decay width 
is effectively suppressed by AqcD/ m b- It is now straightforward to obtain: 



ST 



cqg 



64gjG 2 F \V cb V*\ 2 m b - / dT cqg {2ir)H\m b v - Pl 



■{(<* + C2) [l(fi /Hb +f2 /Hb ) 

+2c lC2 (f^ Hb +ff Hb )}. 



+ 



16 



/ f u/H b f u/H b 
{J 3 + J 4 



P2 ~ k)pi ■ p 2 

i 



(15) 



The parameters W^lp, • • • are defined as: 



fi 



u/H b 



= s 



(u,l) 



(u,l) 
b ' 



f u/H b 
J2 



V 



f u/H b _ <,(ii,8) , ,r,(u,8) f u/H b _ , ,(m,8) a{u,8) 
J 3 - H h + 1 H b > /4 - V H h + A-h , 



(u,l) 
H b 

i,t 
H b 



+ A 



(u,l) 



'3 ~ ^H b H b 

Performing the phase-space integration we obtain: 

<5r cgs = -327ra s r |V;* g | 2 {l + 3z 2 + 2z 3 - 6z 2 ln(z)} 

■{^(c? + 4U /Hb + / 2 U/H 1 + [ 8 -(cj + c 2 ) + 2c lC2 )[/ 3 ^ + fs /Hb }}, 



(16) 



where z = ^% and 

m z 



To = 5^I|F cb | 2 . 



(17) 



(18) 



1927T 3 

Similarly, one can work out the contributions of other two processes. To present our results we 
introduce: 



cug 



cug 



ccg 



1 

3 
1 

6 

21 
3 



[2 - 9z + I8z 2 - llz 3 + 6z 3 ln(z)], 

[11 - 54z + 36z 2 - 2z 3 + 9z 4 - 12z 2 (3 + z) ln(z)], 

l-7z + 6z 2 ] - 2(5 + 2z 3 ) In [l - 3z - (3(1 - z) 
4z 2 (l -/?) 



+4z 3 ln[z(l + /3)] + 10 In 



11 - 86z 



6z 2 + 108z 3 



(1 + /3) 2 
- 10 In 



+ 2 



-9z 4 In 



4z 2 (l-/3) 
(1 + /3) 2 

1 - 3z - 0(1 - z)] + 2z 2 (3 - 2z + 9z 2 ) In [z(l + /?) 



3z 2 + 2z 3 



(19) 
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with j3 = y/1 — 4z. The processes H b (P) — ► c(pi) + u(y»2) + G(/c) + X gives the contribution: 



5T cug = 32na s r \V:f{U[f q 1 /Hb F cug + f^ /Hb F ( 



cug 



+ 



r 4 .2 , 1/2 



4 + ^(c2 + 2cic 2 



fl/Hb TP 

J 3 r < 



, f q/H b p 

cug ~ J4 - 1 cng 



(20) 



Replacing the u quark with a c quark, we obtain the contribution 5r c5s from the process H b (P) —* 
cfri) + c(p 2 ) + G(fc)+X. 
We use the parameters: 

m b = 4.8GeV, « = 0.085. (21) 

Correspondingly we have a s = 0.18, c\ = 1.105 and c 2 = —0.245. The lifetime ratio can be predicted by 
using experimental values as: 



r(H h 



T(B°) 



(B d ) 1 r(H b ) 1 [r(H b )) exp r exp (B°)' 



'T(B )-T(H b 



(22) 



theory 



For our numerical estimation we only take effects of valence quarks into account and use isospin symmetry. 
We obtain: 



r(B-) 
t(B°) 

r(B 8 ) 
t(B°) 

r(At) 
r(B°) 



l + 5.06/i" /B + 7.82/ 2 u /B +8.55/ 3 u/B + 13.26/ 4 U /B ~ 
1 - 0.94/ 1 s/Bs + 1.50/" /B ~ - 2.55/ 2 /Ba + 4.07/ 2 " /B ~ - 



1.59/. 



+2.56/ 3 " /S " 



4.34/ 4 /Bs + 6.94/ 4 u/S " 



1 + 1.36/? /A * + 1.20/" 75 " 



0.69/ 2 /Afc + 3.25/ 2 " /S + 2.27/3 d/Ab 



+2.04/. 



1.23/f Ai + 5.53/£ /B ~ 



(23) 



At moment no information for these nonperturbative parameters is available. This prevents us to give 
numerical predictions. If we use the approximation of vacuum saturation, one may estimate the correction 
to lifetimes of beauty mesons. We use this approximation below to give some numerical predictions for 
mesons. 

If one uses the vacuum saturation approximation, the nonperturbative parameters are either zero or 
can be expressed with the light-cone wave function 0+ of i?-meson, which has been studied extensively(See 
e -g-) HU HH H31 HI] ■ The light cone wave function is defined in the light-cone gauge as [21]: 



(0\q(zlh ■ lTh(0)\B(v)) = --F B Tr [ 757 • ITP V ] <j> + (z), (24) 

where Fb is the decay constant in HQET. At tree-level it is related to the decay constant fs of full QCD 
via Fb = /bV^b- T is an arbitrary Dirac matrix. Under the approximation of vacuum saturation, one 
can show with Eq.(24) that all nonperturbative parameters in Eq.(13) are zero except V^' 1 ^. Hence only 

Q / B 

f? is not zero and it is given by: 



f g/B = f B M B _ J_ J_ 

1 Am b A|' X B 



(25) 
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where (j>+(k) is the Fourier transformed function of <^>+(z). The parameter Xb is at order of \qcd- It is 
estimated to be in the range 0.35 ~ 0,6MeV 20l 124 [ 125 } 126 j . With this we obtain the following numerical 
results in the approximation of vacuum saturation: 



r(B-) 
t(B.) 



1 + (0.15 ~ 0.45) 



1 + (0.017 ~ 0.049) 



(26) 



r(B°) 



where we have used SU{3) isospin symmetry. Comparing experimental results, the predicted ratio of B~ 
is too large. It should be noted that the approximation of vacuum saturation is not a well-established 
approximation, it serves only for a rough estimation. A detailed study of these nonperturbative param- 
eters is required. If one assumes that the nonperturbative parameters of A& is at the same order of 
determined as above, one can also have a large correction to the ratio of A&. 

To summarize: In this letter we have studied corrections of spectator quarks to lifetime ratios of 
beauty hadrons. Formally, these corrections are at order of m^ s . We find some decay channels which 
lead to that these corrections are proportional to certain averages of the squared inverse of the momentum 
carried by a spectator quark. Hence these corrections are effectively at order of m^ 1 . The corrections 
are calculated at tree-level and nonperturbative effects are parameterized with nonlocal operators. Since 
the nonperturbative parameters are unknown, we are unable to give numerical predictions in detail. 
With a simple model one can estimate these parameters and the obtained ratio of B~ is too large in 
comparison with experiment. However, estimations of these parameters in order of magnitude indicate 
that these corrections may be large enough to accommodate experimental results. A detail study of these 
nonperturbative parameters is required to have detailed predictions. 
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